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We construct cosmological models with two scalar fields, which has the structure as in the ghost 
condensation model or fc-essence model. The models can describe the stable phantom crossing, 
which should be contrasted with one scalar tensor models, where the infinite instability occurs at 
the crossing the phantom divide. We give a general formulation of the reconstruction in terms of 
the e-foldings TV by including the matter although in the previous two scalar models, which are 
extensions of the scalar tensor model, it was difficult to give a formulation of the reconstruction 
when we include matters. In the formulation of the reconstruction, we start with a model with 

i^ ", some arbitrary functions, and find the functions which generates the history in the expansion of the 

universe. We also give general arguments for the stabilities of the models and the reconstructed 

("■^ ' solution. The viability of a model is also investigated by comparing the observational data. 

(N' 

, PACS numbers: 95.36.+X, 98.80.Cq 

I. INTRODUCTION 

The observation of the cosmic microwave background radiation (CMBR) shows that the present universe is spatially 
flat [1, 2]. Hence, the starting point for current cosmology is spatially flat Friedman-Lemaitre- Robertson- Walker 
(FLRW) universe, whose metric is given by 

ji^; ds^ = -dt^ + a{tf Yl {dx'f. (1) 

D ■ i=l,2,3 

(—! ■ 
'■^ ■ Here a is called the scale factor. The FLRW equations in the Einstein gravity coupled with perfect fluid are well-known 

to be: 



(N 



>. p = ^H'\ p^^^(3H^ + 2H] . (2) 



X 



'^ I Here the Hubble rate H is defined hy H = d/a and dot denotes the derivative with respect to the cosmic time t. 
^ ■ By the observation of the type la supcrnovae, we now believe that the expansion of the present universe is accel- 

erating [3-5] . What generates the accelerating expansion of the universe is called dark energy (for recent review, see 
CO . [6-8]). The dark energy could be a perfect fluid with negative equation of state (EoS) parameter . The future of the 
^— ^ ' universe is mainly governed by the equation of state parameter wue of the dark energy wbe = Pde/pde, where pde 
^_J , and pDE are the pressure and the energy density of the dark energy. When wde < —1/3, the accelerating expansion 
can be generated and the observational data indicates that wde is close to —1. If the iude is exactly —1, the present 
universe is described by the ACDM model, where the accelerating expansion is generated by the cosmological term 
and the universe evolves into the asymptotic de Sitter space-time. The dark energy with — 1 < wde < —1/3 is called 
as quintessence and that with ujde < — 1 as phantom [9]. If the dark energy is the phantom, the future universe 
Cd . usually evolves to a finite-time future singularity called Big Rip [10], where the scale factor a of the universe will 
diverge in the finite future (for other type of singularities, see [11, 12] and for the classification of the singularities, 
see [13]). 

Recent cosmological data seems to indicate that there occurred the crossing of the phantom divide line in the near 
past (see, for example, [14]), that is, the equation of state (EoS) parameter wue crossed the line wbe = — 1- 

In this paper, we consider the models with two scalar fields to describe the phantom crossing. It is known that 
in one scalar tensor models, the large instability occurs when crossing the phantom divide (or cosmological constant 
border) [15]. We now develop phantom cosmology described by two-scalar tensor theory which represents a kind of 
quintom model [16-18] (for review, see [19] and for generalizations, see [20, 21]). In case of two scalar tensor, we can 
construct a model which is stable at the phantom crossing. In the previous work, it has not been succeeded to a 
general formulation of the reconstruction of the two scalar tensor model when wc include matters. In the formulation 
of the reconstruction, we start with a model with some arbitrary functions, and find the functions which generates 
the history in the expansion of the universe, which complies with observational data (for the reconstruction in F[R) 
gravity, see [22] and the reconstruction for the general models, see [23]). In this paper, in terms of the e-foldings N, 
we consider the two scalar models whose structures are similar to the fc-essence models [24-26] or ghost condensation 
models [27, 28]. In the models, it is straightforward to include matters and it makes us to be able to consider the 
models with the phantom crossing and the matter dominant era. 



II. COSMOLOGICAL RECONSTRUCTION BY ONE SCALAR MODEL AND STABILITY 

Before going to the two scalar model, we clarify the problem in the one scalar model, that is, there appears infinite 
instability at the phantom crossing [15]. This section is based on [21, 29] and we start with the following action: 

d^xV^l^^R-^io{c^)d^cpd'^<f>--V{cP) + cA . (3) 

Here, uj{(j)) and V{(j)) are functions of the scalar field (p. The function uj{(j)) is not relevant and can be absorbed into 
the redefinition of the scalar field cj) as follows. 



'^■ 



d(b^]^W\- (4) 



Then the kinetic term of the scalar field in the action (3) has the following form: 

rj,\a j,oMj, f -dnipd'^^ when w(0) > . . 

-u;{^)d,<pd'^<P^^g^^g,^ whenc.(0)<O ' (^) 

The case of w(0) > corresponds to the quintessence or non-phantom scalar field, but the case of a;(0) < corresponds 
to the phantom scalar. Although U!{4>) can be absorbed into the redefinition of the scalar field, we keep uj{(j)) since 
the transition between the quintessence and the phantom can be described by the change of the sign of oj{<f>)- 
For the action (3), the energy density and the pressure of the scalar field are given as follows: 

p - i^(</))02 _^ y(</)) , p = 10.(0)02 _ y(^) . (6) 

If the potential V is positive, the EoS parameter w is greater than —1 if a;(0) is positive but w is less than —1 if uj{(j)) 
is negative. Then the transition between the quintessence and the phantom can be described by the change of the 
sign of Lo{4>). By using the FLRW equation (2), the EoS parameter w = p/ p can be expressed in the following form 

- ^* (7) 



3ff2 



Then in the quintessence phase, where 1 > l/Zw > —1, we find H < and in the phantom phase, where w < —1, 
H > 0. Then on the point of the transition between the quintessence and the phantom, H vanishes. 

In order to consider and explain the cosmological reconstruction in terms of one scalar model, we rewrite the FLRW 
equation (2) with the expressions (6) as follows: 



2 .•_ 1 



w(0)02 = _J_H, v{(b) ^ ^ {3H^ + H) . (8) 



Assunung uj{4>) and T^(0) are given by a single function /(0), as follows. 



^(0) = -^fA^) , Vi<j>) = 4 (3/(0)' + fA^)) , (9) 

where , = ^ , we find that the exact solution of the FLRW equations (when we neglect the contribution from the 
matter) has the following form: 

= t, F = /(t). (10) 

It can be confirmed that the equation given by the variation over 

= w(0)0 + iw 0(0)0' + 3i/w(0)0 + VA(b) , (11) 

is also satisfied by the solution (10). Then, the arbitrary universe evolution expressed hy H = f{t) can be realized 
by an appropriate choice of aj(0) and 1^(0). In other words, by defining the particular type of universe evolution, the 
corresponding scalar-Einstein gravity may be found. 



We now show that there occurs the very large instability when crossing the cosmological constant line w ~ —1. For 
this purpose, by introducing the new variables X^f, and Y as follows 

Z^^, Y^lf, (12) 

we rewrite the FLRW equation (2) with (6) and the field equation (11) as 

d^__fMM^^^_.Uy_v^ dV _ Ui^)il-ZY)Z 

dN~ 2f.4cP)H '^^ '' dN^ H^ • ^'^' 

Here N is called as e-foldings, defined by the scale factor a as a = c^^^", where A'o is the value of N in the present 
universe. Then we find d/dN = H~^d/dt. Since we have Z — Y — 1 for the solution (10), we now consider the 
following perturbation: 

Z=1 + SZ, Y = 1 + 6Y (14) 

Then 

Here, the solution (10) is used. The eigenvalues M± of the matrix are given by 

In order that the solution (10) could be stable, all the eigenvalues M± must be negative. One now considers the 
region near the transition between the quintessence phase and the phantom phase, where H ^ 0. When we consider 
the transition from the quintessence phase, where iJ < 0, to the phantom phase, we find H > 0. On the other hand, 
when we consider the transition from the phantom phase, where H > 0, to the quintessence phase, we find H < 0. 
Then for both transitions, one finds H/HH is large and negative. The eigenvalues M± are given by 

A/+ ^, M_-0. (17) 

HH 

Then M+ is positive and diverges on the point of the transition H = 0. Hence, the solution describing the transition 
is always unstable and the instability diverges at the transition point and therefore the transition is prohibited in the 
one scalar model"'^. 

III. COSMOLOGICAL RECONSTRUCTION BY TWO SCALAR MODEL 

We now consider the new type of two scalar models and give a formulation of the cosmological constant. The two 
scalar models are different from those proposed in the previous papers [16-18] and the action with two scalar field (f> 
and X is given by 




S 



J d^xV^ (-^R + X+ uj{4>)X^ -U + 7]{x)U^ - V{cl), x) + C„?j , (18) 



and assume the flat FLRW metric (1). Here X and U are defined by X = —^d'^(f)d^(j) and U ~ ^^d^^xdfiX and Cr, 
is the Lagrangian density of the matters. 



^ This does not always mean that any one scalar model with non-canonical kinetic term shows the instability (see [30], for example). 



We consider the models where there do not exist no direct interactions between the scalar fields and the ordinary 
matters. Then the conservation law of the matters is given by 

p;,(7V) +3(p™(7V) +p„AN)) =0. (19) 

Here p,„ and Pm are the energy density and the pressure of the matters, respectively and the prime ' denotes the 
derivative with respect to the e- folding number N . 

For the action (18), the FLRW equations have the following form: 

X + 2uji^)X^-U + 2rj{x)U^ = -^HH' + \p'^. (20) 

The energy density contributing to the Hubble rate can be divided into contributions from the scalar fields and the 
matter part as follows, 

i?'(iV) = y(p.(iV)+Prn(iV)), (21) 
and we rewrite the FLRW equations as follows, 

F(0, x)~\{X-U)^ p,{N) + ip;(iV) , (22) 

W(0)X2 + 77(x)t/2 + \{X -U) = -^Ps{N) . (23) 
The field equations for the scalar fields are given by 

1/0(0, x)+3tJ,0(</-)x2 ^ -H'^{l + QujX)(t)" - {HH'{1 + &ujX) + ZH'^ + QH'^LoX)(j)' , (24) 

V^{<P, x) + HxixW^ = H'-il - 677[/)x" + {HH'il - QijU) + 3H^ - (SH\U) x! , (25) 

where </> = ^ and % = ^ respectively. Then, if V, to, and 77 satisfy the following relations 

V{mN, mN) = /(TV) + -^f\N) , 
Uj{mN) + riimN) = -- ■' ^ ' 



im'^H^iN) ' 



V.^{mN, mN) + -m'^H^{N)LU,4,{mN) 



= -m{3H^{N) + HH'{N)) - 3m^H^{N)u{mN){H^{N) + HH'{N)) , 
3 
4' 
m{3H^{N) + HH'{N)) - 3m^H^{N)T]{mN){H^{N) + HH'{N)) , (26) 



l/x(miV, mN) + -m^H^{N)ri^^{mN) 



we can obtain the following solution 



X = miV, H\N) = —{f{N)+p^{N)). (27) 



We now find explicit forms of the functions a;(0), rj{x)i and V{(j), x) in the action (18), which satisfy the above 
relations. We may choose, for example, the forms of w((/)) and rj^x) with two arbitrary functions a(0orx) and 
/(0/77iorx/?Ti) as follows, 

Vix) = -a(x)- (29) 

We also define a new function /(</>, x) by 

J{<P,X) = -m[ [ d<p{Zm^uj{(l>)H\<Plm) + 2H^{^/m)) 
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+ Jdx {3m'vix)H^ ix/m) - 2H' (x/m)) | , (30) 

which gives 

fimN, niN) = f{N) . (31) 

Using the function /(0, %), we find V{(l), x) is given by 

l/(0,x) = -m{ I d4,H''{4>lm)- fdxHHxM 



By reconstructing the action by using two functions, a(0orx) and f(<j)/movx/m)^ wc can obtain the arbitrary history 
of the expansion of the universe as wc desire. 

IV. STABILITY 

In this section, we discuss the stabihty of the reconstructed solution by considering the perturbations. We also 
investigate the stability by using the sound speed. If the square of the sound speed is negative, the fluctuations grow 
up exponentially and the homogeneous universe becomes unstable. 

A. Stability of the reconstructed solution 

In order to investigate the stability of the solution (27), we consider the fluctuation from the solutions. We define 
the 1st order perturbations from the solutions as follows, 

== 00 + <^0(iV) , x = Xo + 5x{N), = 0o + 5x(7V), x^Xo + 5v{N), (33) 

where 5x{N) = 5(j){N) and 5y[N) = 5x{N). We note that this system has four physical degrees of freedom (dof) 
and we identify 0, X: </*: X with the dof here. It is also noted that the Hubble rate H are not the dof because it can 
be removed by the constraint equation, i.e., the first FLRW equation (21). Substituting the expressions in (33) into 
Eqs. (22), (23), (24), and (25), we obtain the evolution equations of the 1st order fluctuations as follows. 





Here 



if-i 
H-^ 

ai 02 03 04 



_2HH' - 2H'^ + Sm^uj^^H^H' + m^K^H^{l + m^uH^) 

+ 1 (aHH' + 3if'2 + HH" + !^(1 + m'ujH^){H^ + HH') 
H \ 2 



(34) 



"^ - i/(lV3^W) W"" + ^^')(-^ + ^"^^^ ) - ^ ) ' 



(1 + m^UjH^) ^ ;r^ 



03 = - 3 (1 + rn^LoH^) 






2 1 + Sm^ujH^ 



H{-l + 3mhjH^) \2 



-{H^ + HH'){1 + Sm^ujH^) + H' 



ae = 



2HH' - 2if '2 - im^ii^H^H' + m^K^H^{-l + m^riH^) 
H{-l + 3m^T]H^) 

- (aHH' + 3H'^ + HH" + !^(-l + m'7^H^){H^ 
H \ 2 



HH') , 



a? 



m'^K'^ —1 + ■m'^riH^ 



-1 + 3m'^i]H^ 



(l + Sm^wff^), 



as = - 3 - 



-{-1 + m^riH^)- 



6m''f]H^H' + 'im^^Ji^H^ 



2 ' ' '" ' i7(-l + 3m2?7iJ2) 

where we have used the constraint equation between the perturbations, which is obtained from (21) as follows, 



(35) 



5H = 



6H 



= 0(1 + 30.02)^ + ( ^C.,^^4 ^ y\ g^ ^ ^(_^ ^ 3^^2)^y ^ (|^,xx' + ^x ) -^X . 



(36) 



In order the solution (27) is stable, the real parts of all the eigenvalues of the above matrix (34) must be negative. 
Then we find the conditions that the real parts of all the eigenvalues are negative. The characteristic equation and 
its solutions are given as follows: 



A^ +AiX^+A2\''+A3X + A4 



' =2 



0, 
Ai 1 



-9±v/92-4S±4pj-^, ^(-g±Vg2-4S±4p)-^. 



(37) 



Here 



Ai ~ -as — as, A2 = 0303 - 0407 — 



ai + og 
H 



1 r 



asttQ - 0405 - 0207) , A4 = — 2(0106 - 0205) . 



Here S is one of the solution of the following cubic equation; 



"3 



2 



S2-S3S 



^ , B1B3 



-Bl^Q, 



that is 



Here 



S= -P-v/P' + Q- 



Si 
6 



Si 



iAl 



8 
53 



^2 



B. 



AI A1A2 



^3, 
- - 8 2 

p ^ ^1 . B1B3 Bl o = ^^ ^^ 
216 6 16 ' 36 3 



-P + v/p2 + 03 ' 

3At , ^2^2 A1A3 



(38) 



(39) 



(40) 



B.3 



256 



16 



A4 



In (37), (p, q) are defined by 



VS - B3 , 



v/2S - Bi 



(41) 



(42) 



Now we consider the real solution of the cubic equation (39) since we like to find the condition for the maximum value 
of the real parts of A (denoted as SRAmax) to be negative. Then, all the quantities given by the capital characters in 
(38) and (41) become real numbers. Then we find the expression of S for four cases. 



2|g|2cos 






6 



^ = < 



P+^/WTQ' 



p-^/W 



for g < and P^ < \Q^\ 
^ for g < 0, P2 > |q3| gjj^j p > 



-P + /P^ 



-P-^W+Q- 



Bi 
6 



for g < 0, P2 > |g3| gjj^j p <- 



- p+ /p^Tg^ 



-p + v/s2 + g3 ■ + ^ for g > . 



We also find the expression of SRAn 
For S^ > Bs, 



5RA„ 



i (g + V-2S -Bi+ 4p) - 4^ for S > ^ and - 2S - Bi + 4^ > 

ig - 4^ < for S > ^ and - 2S - Bi + 4? < 

^y/-2E-Bi+Ap - 4^ for S < f' and - 2S - Si + 4/ > 

-4^ forS< 4^ and -2S-Bi+4/ <0 



<0, 



for S2 < B3, 



■JT-^max ~^ 



Ul+\ 1 ( -22 - Si + V(25 + Si)' + 16p2 ) - 41 forS> 



2 ^ i (-2S - Bi + v/(2S + Bi)2 + 16p2 j _ 4^ forS<f 

< . (43) 

This condition may reduce to the constraints on a' [mN), a(mN), f{N) in ui^^imN) etc. However, it is not straight- 
forward to obtain the constraints on a'{77iN) for each of models (or equivalently, for each of f{N)) explicitly because 
the condition, P^ > \Q^\, is the 12th order inequality with respect to a'{mN) and we cannot solve it analytically in 
general. Furthermore, as in the case of the one scalar k-essence model [31], it is very difficult to find the explicit form 
of a{mN) even if we obtain the explicit constraints of a'{mN) since it is not an equality but an inequality. Therefore, 
from the condition for the stability of the solution, it is difficult to find the explicit form of the function a{mN) for 
each of models. We can obtain, however, more explicit constraint for a{mN) from the condition that the square of 
the sound speed should be positive, as we will see in the next subsection. 

B. Constraints from sound speed 

We now consider the stability constraints coming from the sound speed of the scalar fields. By considering the 
perturbation of the scalar fields, we find the expressions of the sound speed c^^ (« = (^, x) &s follows. 



„2 _ P't>,x _ 1 + 2u)X 2 _ Px-U _ -1 + 27?t/ 
^ " l + QuX ' ""-'^ " J~^ " -1 + 6?7C/ 



^?. = ^ = ^^l^. 4 = ^ = ^^^f^- (44) 



The sound speed must be < c^j(< 1) for the stability. This condition can be expressed as the constraints for the 
function a for the reconstructed solution (27), 

am,=mN>^^, or a(^)|,_^<_^-_l-^, 
and a(x)lx=mjv > or a(x)lx=mw < TEH- (45) 

Although the models in this paper could be regarded as an effective and classical low energy theory of more 
fundamental theories, we may also consider the (in)stability as a quantum theory. As a quantum theory, we require 
the conditions P4i,x > 0, /O0,x > 0, p^^u > 0, and p^^jj > for stability and the conditions give the constraints 
a{mN) < — ^'2jj-2 and a{mN) > 3^4 ^4 ■ If the universe has the phantom crossing time, however, the derivative of the 
energy density for the scalar fields /' in (27) changes its sign from negative to positive and therefore the reconstructed 
solution (27) of the evolution of the universe cannot be stable as a quantum theory at the phantom crossing time. 
We should note, however, that this analysis of the sound speed is valid only in the sub-horizon scale and in absence 
of the gravity because we do not consider the FLRW space-time and the gravitational fluctuations. Therefore, we do 
not take care of the quantum instability here and we only consider the classical stability of the models in the next 
section. 



TABLE I: The constraints for the cosmological parameters from the cosmological observation [32]. 



observables 


68% CL 


WOEO 


-1.032 ±0.144 


10 


(-0.60, -0.30) 


jo 


(-0.88, 0.90) 


So 


(-0.57, 1.07) 



V. TWO SCALAR MODELS 
A. Model 1 

As a first example of tlic model with phantom crossing, we consider the following model with matters: 

H-2/^^^ _ rr2 /^o C0sh[7(iV - A^^)] , ^ 3(^N-No)\ ^fi^ 

H (N) - H, l^n^^,—^^^-^-——^^ + J7,„oC J , (46) 

noEo =0.74, r!™o = 0.26, (47) 

where subscript "0" denotes the present value of each parameter and the free parameter Nc is the e-folding number 
at the phantom crossing time. The parameter 7 is to be determined by the observational data and the stability 
constraints. The model (46) reduces to ACDM model in the limit of 7 = 0. 
At the point of the phantom crossing N = Nq, we find 

H = H„ H = 0, H= ^I^W^. (48) 

2cosh[7(7Vo-Afe)] ^ ' 

Therefore if we try to realize the development of the universe given by (46) by using one scalar model (3), the 
eigenvalue M+ in (17) diverges at the crossing point and therefore the infinite instability is generated there. This tells 
that the development in (46) cannot be realized by one scalar model. In case of two scalar model, however, as wc will 
see soon, the reconstructed solution becomes stable and the universe can develop as in (46). 

In order to find the viable region of the parameter 7, we introduce cosmological parameters, the equation of state 
parameter wde, the deceleration parameter g, the jerk parameter j, and the snap parameter s, whose present values 
are constrained by observations [32]: 



/' 1,3 kV ^d^a . Irf^a 

6J a 2 6//q a dt-^ a dt^ 



The observational constraints to these parameters are shown in Table I. In case of our model, the observable parameters 
can be expressed as 

WDEO = -1 - I tanh[7(Aro - N,)] , 

3 7 

9o = -1 + 2^™o - ^^DEo tanh[7(Aro - N^)] , 

.70 = 1 + f^DEo Q7 + 37tanh[7(7Vo - N,)] 

So = 2 ("deo7^ tanh[7(Aro - N^)] - 27r!,„o) + (jo + S^o + 2)(3 - 90) + iql - 2 . (50) 

We show the viable region of the parameter 7 for various values of Nc in Table II. In order that the observational 
constraints could be satisfied, the phantom crossing must occur in future for any model (not in past), which is clear 
from the form of parameter q in (49). When we change the variable form the cosmic time t to the e-folding number N, 
we used the Hubble parameter in ACDM model as an approximation since 7 is nearly for each case and therefore the 
energy density of the scalar fields do not vary so fast. Then the relation between the cosmic time t and the c- folding 
number N is now given by 

, ^ 2 , ^ /■*«,, 2 sinh[criol , ^ 

H{t)o^-acoth[at], No - N = I H{t)dt c^ -\n—^ , (51) 



TABLE II: The viable region of the parameter 7 for various values of N^ 



N,-No 


tc — to 


l7l 


1.0 


~25.2Gyrs 


0.236 < 


l7l 


< 0.306 


5.0 


~ 150.2 Gyrs 


0.101 < 


l7l 


< 0.129 


10 


~ 300.7 Gyrs 


0.074 < 


l7l 


< 0.096 



where (7 ~ -^-v/SA. We assume the universe was dominated by the matter from the time t,„. When to — 13.7 x lO^yrs 



and tm = 7x10"* yrs, we obtain Nq — N„ 



.16. If we assume that the expression in (51) is vahd when Nc — Nq < 10, 



we can estimate the crossing time tc for each Nc — Nq as shown in Table II. 

As a second step, we try to find the function a{mN) so that the reconstructed solution (27) could be stable. We 
investigate two cases, a{mN) = — m'^H'^ ^^^ a{mN) = + ^-iff-i by putting m = —. For both of the cases, we find that 



the solution is stable for ten or hundreds billion years (depending on the value of N^) after the phantom crossing but 
becomes unstable at a later time. See Figure 1 and Table III. 



Stability of the solution 




FIG. 1: 3?A„ 



(A^ — Nc) curve in the case of a(mN) 



— jW and 7 — 0.306. The region KA„ 



< is stable but ^Xmax > 



is unstable. This graph shows model 1 is stable for a few e-folding number after phantom crossing when N — Nc = 0. 



TABLE III: The limit values of A — Ac and t — tc when model 1 is stable in the case a{mN) = — ttti- We obtain similar result 
also in the case a{mN) — ttt^- 



Nc - No 


Stable limit 


I7I 


1.0 


A- 


-Nc 


_ I 1.253 
— +1.12 


t-tc^ +i!:i Gyrs 


|_| 0.236 
|7| -0.306 


5.0 


A- 


-Nc 


,5.524 
— +5.402 


i-tc - 


- +119.1 Gyrs 


U,l 0.101 
!7l =0.129 


10 


N- 


-Nc 


1 10.569 
— +10.440 


t-tc^ 


- +117:3 Gyrs 


U,l 0.074 
|7| -0.096 



Finally, we explicitly give the reconstructed Lagrangian of the scalar fields which generates the stable phantom 
crossing in the case m = - and a — 



17W- 



£., =X 



K'^jH^noEQ sinh[K7((/) - (/i^)] k'^A 



cosh[K7(0o - <j)c)] 

Hq i^DEO 

H? \cosh[K-f{(j)Q - 4>c)] 
H^n,no{A + 2] 



H^inq 



3- 



3A + 2 



4k2 



g-3K(0-0o) 



cosh[K7( 



X^-U - 



_ifA_ 

H^Kx) 



U' 



^ ,, ,'37 3A + 



smh[Kj{(j) - (j)c)] 
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TABLE IV: The viable region of the parameter S for the model 2 and e for the model 3 for various values of A'^c 



N,-No 


tc — to 


5 or e 


1.0 


~25.2Gyrs 


0.37 <S OT e< 0.89 


5.0 


~ 150.2Gyrs 


0.14 <6 OT e< 0.28 
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~ 300.7 Gyrs 


0.10 < <5 or e < 0.26 



1 / H^noEa \fiA + 2 3A + 3 

— cosh[K7(x - Xc)\ H sinh[K7(x - Xc) 



2 \ cosh[K7(xo - Xc)] 

3k(x-Xo) 



-ffo^mo(^ + 2 



4k2 



(52) 



where A = 1, 



B. Model 2 

As a second model, we consider the following model which may have the Little Rip [33-40]: 

H [N) - gp i^DEo ^^^_^^ ^ ^^^^ ^ ^^^ + »moe ^J (53) 

where 5 is the parameter which should be constrained. The obtained constraints are given in Table IV. Then, as in 
model 1, we choose a{mN) — — sijlif lyi (where A = 1, —2) and we obtain results similar to those for model 1 about 
the stability, that is, the solution becomes unstable in ten or hundreds billion years after phantom crossing. This 
means that the universe may evolves to other future without the Little Rip singularity for a{mN) — ^y^ij^i/j^) ■ The 

reconstructed Lagrangian of the scalar fields in the case to = i and a = ^i'jji is given as follows: 

{c-'^^^'l'"-'!'^) + k6{(I)o - (t>c) H-^{K(f>) m{n(t))) H^i^x) 

1/ H^noEo \(f, 3A + 2 3A + 3 3<5\ _,„... 3k'S{A + 1)^_^ ^,2 



,.2 



(;-k5(0o-0c) _|_K,5(0o - 0c)/ I V 4 (5 4 

.6 (3 - ^) (0 - 0.) + ^} - M^l.^|^e-(.^.o) 
1 / H^noEo \(f3A + 2 3A + 3\ _,5(v-v ) , 3^25(^+1) 



-''---'^+ '"'"':^'\ x-xcf 



^*^(X - Xc)U ^%l(A±^e-(x-xo) . (54) 



.2 yg-K<5(xo-xc) + K(5(xo - Xc)/ I V 4 S 

4«;2 

C. Model 3 

As a third model, we suggest the model which may have asymptotically de Sitter space. The model has the following 
form: 

e-4N~N^} ,^(l_ -e{N-N^)\ 
H\N) = go^[»z..o ^,,(^„.^4|;^_^_,(^„_^^;^ + a„oe-(---")] (55) 

where e and 9{< e) are the parameters which should be constrained. We set 6 = 0.01 for simplicity. In this case, we 
obtain the same result as in the Little Rip model (53), which is shown in Table IV. As former models, we choose 
a{niN) = ^2^(jv) (where A = 1, —2) and we obtain the result similar to former models (46) and (53) about the 
stability, that is, the solution becomes unstable in ten or hundreds billion years after phantom crossing. This tells that 
the universe may evolves not to asymptotically dc Sitter space for a{mN) = — twttw\- The reconstructed Lagrangian 

of the scalar fields in the case ra = - and a = ,„, is given as follows: 



K^eH^noEo e'i-<'t>-'t>.) _ e-«e(0-0<=) ^2^ \ ^2^ 



'^c-'=<'^o-'^<^) + f(l-e-'=^W«~'^<=)) H^{K(t>) H^K(j))J H^ii^x) 
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1 



'-o'-l'DEa 



1 

K^ \e-'^<^('t' »-'!><:) + 1(1 — Q-K.f.{cj>0-(l>a)\ 

3A + 2 3A + 3 36* 



3 _ 3A + 2 ^ 3A + 3 _ 3e\ ^_.,(^_^^ 
4 e 4 



io>'!'DEQ 



e f /3A + 3 3A + 2 



4 y 4 



^3A + 2 3^ + 3\^_(,_,^) 



V 4 e 

3A + 2 



"j c«^(x-x<=) + (3^ + 3)kx + 



ggr^mo(^ + 2) 3^(0_0p) 

4aj2 



gg»„o(A + 2) 
4k2 



-3k(x-Xo) 



(56) 



VI. SUMMARY 

In this paper, we considered the models with two scalar fields, which has the structure as in the ghost condensation 
model [27, 28] or k-essence model [24-26]. The models can describe the stable phantom crossing, which should be 
contrasted with one scalar tensor models, where the large instability occurs at the crossing the phantom divide. In 
the previous two scalar models, which are extensions of the one scalar tensor model, it was difficult to give a general 
formulation of the reconstruction when we include matters, which is realized in this paper in terms of the e-foldings N. 
We also gave general arguments for the stabilities of the models and the reconstructed solution. We also investigate 
the viability of the models by comparing the observational data. 
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